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Abstract 

We discuss a possible answer to the fundamental question of why nature would actually 
prefer low-scale supersymmetry, but end up with a supersymmetry scale that is not com¬ 
pletely natural. This question is inevitable if we postulate that low-energy supersymmetry is 
indeed realized in nature, despite the null observation of superparticles below a TeV at the 
Large Hadron Collider. As we argue in this paper, superparticles masses in the multi-TeV 
range can, in fact, be reconciled with the concept of naturalness by means of a cosmological 
selection effect—a selection effect based on the assumption of an exact discrete /7,-symmetry 
that is spontaneously broken by gaugino condensation in a pure supersymmetric Yang-Mills 
theory. In such theories, the dynamical scale of the Yang-Mills gauge interactions is required 
to be higher than the inflationary Hubble scale, in order to avoid the formation of domain 
walls. This results in a lower limit on the superparticle masses and leads us to conclude 
that, according to the idea of naturalness, the most probable range of superparticle masses 
is potentially located at the multi-TeV, if the inflationary Hubble rate is of O(10 14 ) GeV. 
Our argument can be partially tested by future measurements of the tensor fraction in the 
Cosmic Microwave Background fluctuations. 



1 Introduction 


The non-discovery of any supersymmetric partners of the standard model particles (sparticles) at 
the Large Hadron Collider (LHC) experiments so far has excluded sparticle masses in the range 
of a few hundred GeV [Ij. Besides, the observed Higgs boson mass of a 125 GeV [2] suggests that 
sparticle masses most probably lie in the multi-TeV range [3J (see also |3j for a recent analysis), 
in case supersymmetry (SUSY) is really realized in nature. For such a high-scale SUSY scenario, 
serious questions, however, arise regarding the qualification of SUSY as a solution to the fine- 
tuning problem in the Higgs potential. Why would nature prefer low-scale SUSY, but end up 
with a SUSY scale that is not completely natural? What are the physical constraints preventing 
the SUSY-breaking scale from being lower and, hence, perfectly natural? 

In this letter, we show that these fundamental questions could be possibly answered if the 
scale of the last inflation is very high. As we are going to argue, the key element to a better 
understanding of a high SUSY scale is the domain wall problem related to the spontaneous 
breaking of a discrete ^-symmetry in the early universe. Since the formation of domain walls 
after the end of inflation is disastrous for the habitability of the universe, any given inflation scale 
implies a lower bound on the scale of i?-symmetry breaking [5]. Meanwhile, the SUSY-breaking 
scale (and hence the masses of sparticles) and the i?-breaking scale are strictly tied to each other 
by virtue of the flatness condition of the universe. As a result, invoking cosmological selection 
for habitable universes, we find that the probable range of sparticle masses deduced from the 
idea of naturalness can indeed lie at the multi-TeV, if the Hubble rate is of O(10 14 ) GeV. 

The organization of this paper is as follows. In section [2j we discuss the constraint on the 
i?-breaking scale resulting from the cosmological domain wall problem. In section [3j we then 
discuss the probable range of sparticle masses according to the concept of naturalness. The final 
section is devoted to summary and conclusions. 

2 Spontaneous Breaking of Discrete R-symmetry 

For low-scale SUSY to be realized in nature, the SUSY-breaking scale, | 1//2 , should not be 
too large. In addition, the vacuum expectation value (VEV) of the superpotential, Wq, is also 
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required to be small, so as to achieve an almost vanishing cosmological constant, i.e., 


Vo = \F\ 2 - -i- |W 0 | 2 ~ 0 , (1) 

iw PL 

where Mpl denotes the reduced Planck scale. Thus, in order to realize low-scale SUSY, it is 
inevitable to invoke a symmetry—an ^-symmetry—which prevents Wo from being very large. 
A small, but non-vanishing VEV of the superpotential can then be provided by the spontaneous 
breaking of this i?-symmetry0 

Let us emphasize here that global symmetries are generically believed to be violated by 
quantum gravity. Therefore, if we require that our A-symmetry is an exact symmetry, it must 
be a remnant of a gauge symmetry, which appears as a discrete symmetry in the low-energy 
effective theory. In the following argument, we will have a very strict attitude towards global 
symmetries, rejecting them as possible candidate symmetries suitable to protect the VEV of 
the superpotential. Instead, we will assume that the /?,-symmetry protecting the VEV of the 
superpotential is an exact discrete symmetry, Znr. 

A serious drawback when invoking an exact discrete symmetry and its spontaneous breaking 
is the domain wall problem [6[ [7j. If the Hubble rate during the last inflation, H , is very high, 
spontaneous R breaking takes place after the end of inflation, which is accompanied by the 
formation of domain walls. Since we assume an exact discrete symmetry, these domain walls 
are stable and they immediately dominate the energy density of the universe once they are 
formed. Therefore, we need to require that i?-symmetry breaking takes place before/during the 
last inflation, which leads to a constraint on the scale of iLsymmetry breaking. 

To see how the if-symmetry breaking scale is constrained, let us consider a pure SUSY Yang- 
Mills theory, which is, in fact, the simplest model exhibiting spontaneous ii-symmetry breaking. 
In pure SUSY Yang-Mills theories, i?-symmetry is spontaneously broken by dynamical gaugino 
condensation [8j. For example, in an SU(N C ) Yang-Mills theory, discrete Z 2 n c r symmetry is 
spontaneously broken down to Z 2 R symmetry and the resulting Wo is roughly given by 

W 0 ^ A 3 , (2) 

with A denoting the dynamical scale and where we have omitted allowed numerical coefficients. 
In the early universe, spontaneous i?-synunetry breaking takes place when the Hubble parameter 

L It is logically possible that a small Wo is obtained by fine-tuning. Such a possibility is, however, retrogressive 
to the naturalness argument. See our discussion at the end of section[3l 
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H (or the temperature of the universe T) drops below the dynamical scale, H A (or A). 
Therefore, to avoid the formation of domain walls after inflation, the dynamical scale is required 
to be higher than the Hubble rate during inflation, i.e., H [5j, which puts a lower limit on 
the scale of 77-symmetry breaking. 

In Fig.JTJ we show the corresponding lower bound on the gravitino mass, where 717 . 3/2 = 
Vto/Mp L , for a given Hubble parameter during inflation, 

™-3/2 > m 3/2( H ) = - 10 °) TeV x ( 10 i4G e v ) ' ( 3 ) 

Here, the range 0(1 — 100) reflects our ignorance of the precise relation between Wo and A as 
well as between the critical time during the 77-symmetry breaking phase transition and A. To 
obtain a more precise constraint, further quantitative understanding of the non-perturbative 
aspects of SUSY Yang-Mills theory is necessary]^] To give a representative example, we show the 
constraint by taking the relation Wo — A 3 literally, although we should keep in mind the above 
uncertainties. In the shaded region, the corresponding dynamical scale is large enough so that 
17-symmetry breaking takes place before the end of inflation. Therefore, there is no domain wall 
problem in the shaded region. 

According to the above argument, we find that 7773/2 below the TeV range is prohibited in con¬ 
sequence of the domain wall problem, unless the Hubble scale is much lower than 0(1O 14 ) GeV. 
Interestingly enough, such a large Hubble rate corresponds to a rather large fraction of tensor 
modes in the primordial fluctuations seen in Cosmic Microwave Background (CMB), 

/ TJ \ 2 

r ~ 0.16 x [ ——I . (4) 

\10 14 GeV ) v ’ 

Future experiments such as CMBPol |12| and LiteBIRD [13] are expected to reach values of the 
tensor-to-scalar ratio r of 0(1O -3 ). This, thus, opens up the possibility to put a stringent lower 
limit on the gravitino mass through CMB observations, assuming that the spontaneous breaking 
of 77-synunetry is accounted for by gaugino condensation. 

Before closing this section, let us remind ourselves that there are possibilities to evade the 
domain wall problem, even if we assume an exact discrete 77-symmetry. For example, let us 

2 If we rely, for example, on so-called naive dimensional analysis ng the right-hand side of Eq. |2j) is sup¬ 
pressed by a factor of 0{( 47t) 2 ). In that case, the corresponding lower limit gets weaker by the same factor if the 
condition is given by A > H. 
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Figure 1: Constraint on the gravitino mass for a given Hubble rate during inflation. We assume that 
R-symmetry breaking is caused by gaugino condensation as in Eq. |2|). In the shaded region, the cor¬ 
responding dynamical scale is large enough so that R-symmetry breaking takes place before the end of 
inflation. We show the boundary of the shaded region by a thick line to warn against our ignorance of 
the exact relation between Wq and the critical time associated with the R-breaking phase transition. 


consider a model of ^nR-symmetry breaking where a gauge singlet 4> with R-charge 2 possesses 
a superpotential of the following form, 

\n 

W = v n 4> +-< b n+l + ••• . (5) 

n + 1 

Here, v n and A are parameters and we take the reduced Planck mass to be unityj^] The ellipsis 
denotes higher-dimensional terms in (j) which are consistent with the Z 2 n R symmetry. In the 
vacuum, .R-symmetry is broken by the VEV of q i, leading to the VEV of the superpotential, 

Wo - v n+1 /\ . (6) 

In this example, it is always possible to avoid the domain wall production if the singlet obtains 
a large negative Hubble mass squared. Including a negative Hubble mass squared, R-symmetry 
is forced to be broken during inflation. In this case, there is no constraint on the gravitino mass 
from the requirement of no domain wall formation after inflation, since R-symmetry has already 
been broken during inflation^] 

4 The parameter v n should be small to explain the smallness of Wo- Such a small v n could be achieved by 
assuming some additional dynamics behind v n . 

4 One might wonder that the domain wall might be formed once <j) starts moving around its origin after inflation. 
However, the formation does not take place for n > 3 El. 
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3 Naturalness and Sparticle Masses 


Let us now discuss how the above observation enables us to answer a fundamental question 
brought upon us by the results of the first run of the LHC: why would nature prefer low-scale 
SUSY, but end up with a SUSY scale that is not completely natural? For that purpose, let us 
first review the conventional argument on the “natural” range of the SUSY breaking scale and, 
hence, sparticle masses, msuSY- To discuss natural ranges of parameters, it is often transparent 
to consider an ensemble of vacua (or theories) with various SUSY-breaking scales. Here, we 
call the ensemble of vacua the landscape of vacua , adopting the terminology coined in Ref. m 
having a string theory landscape in our mind [15] (see also m for an earlier discussion). One 
way to find the range of msusY preferred by the concept of electroweak naturalness is to consider 
the distribution of different values of the electroweak symmetry breaking scale, uew> f° r a gi yen 
msusY- That is, we should collect all vacua with the same value of msusY from the landscape and 
count how often we respectively encounter each value of uew- Then, for a given value of msusY, 
the distribution of uew is expected to peak around uew — t3USY, since electroweak symmetry 
breaking is triggered by the sparticle massesj^] Thus, given the observed electroweak scale, 
uew — 174 GeV, we infer that, from the standpoint of electroweak naturalness, the sparticle 
masses are most likely to lie in the range of a few hundred GeV . 

An alternative way to deduce the range of msusY from electroweak naturalness is to consider 
the distribution of msusY for a fixed value of the electroweak scale, instead. To illustrate the 
idea behind this alternative, let us start from the initial distribution of msuSY in the landscape, 
without imposing any cuts on the ensemble of vacua ((a) in Fig. [2]). Here, we assume that the 
prior distribution is not severely biased towards large values of msusY, although we do not 
need to know the exact distribution. Since we are interested in vacua with an almost vanishing 
cosmological constant, we restrict the landscape in the next step to vacua corresponding to 
an (almost) flat universe. In this restricted landscape, we expect that the distribution is now 
sharply biased towards low energies, since a flat universe can be achieved more easily for lower 
values of the SUSY-breaking scale [17] ((b) in Fig.[2]). Finally, we restrict the landscape further, 
so that all vacua in the landscape have the same electroweak scale uew- Then, the distribution 
of sparticle masses is cut off around msusY — «EW) since electroweak symmetry breaking is 

throughout this paper, we assume the Minimal Supersymmetric Standard Model (MSSM), so that the Higgs 
potential is given by known couplings and soft breaking parameters of O(msuSY). 
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(a) Landscape of vacua 


(b) Landscape of vacua 
of the flat universe 


(c) Landscape of vacua 
of vew — 174Ge V 



VEw~174GeV 


Distribution of m S usY 


msusY 



Figure 2: Schematic illustration of how to deduce the probable range of msusY- In each panel, the dots 
denote vacua. From (a) to (b), we restrict the landscape to vacua with an almost vanishing cosmological 
constant, which, as we expect, leads to a sharply peaked distribution of m susy biased towards low 
energies. From (b) to (c), we restrict the landscape to vacua with an electroweak scale of Vew — 174 GeV, 
which results in a distribution peaked around vew- 


triggered by the sparticle masses in the MSSM ((c) in Fig. [ 2 ]). As a result, we end up with a 
distribution of msuSY which peaks around uew- This means once again that sparticle masses 
most probably he in the range of a few hundred GeV in view of electroweak naturalness. 

In both approaches, we end up with more or less the same conclusion: msusY should be at 
around a few hundred GeV. (In the following, we shall call the former approach the frequentist 
approach and the later approach the Bayesian approach.) Here, it should be emphasized that 
our deductions of the range of msusY are based on the principle of mediocrity [18], i.e., in both 
approaches we have assumed that we are typical observers living in a typical habitable vacuum. 
This is the reason why we obtained a similar conclusions in these two different approaches. 
The big problem now is that the resultant ranges of msusY in both approaches are in tension 
with the null discovery of sparticles at the LHC as well as with the observed Higgs boson mass. 
This means that, unless we find a way to depart from the above argument and alter the ranges 
deduced above, we almost lose ground on the postulation of low-scale SUSY from the standpoint 
of electroweak naturalness. 

The above conclusion should, however, change if there are crucial restrictions not accounted 
for in the above argument. In the following, we are going to argue that the domain wall problem 
related to i?-symmetry breaking corresponds exactly to such a missing selection rule. Since 
we are trying to shave the distribution function of msusY by applying the missing selection 
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rule for habitability along with other habitable conditions, the cosmological constants and the 
electroweak breaking scale, it is more transparent to take the Bayesian approach. One caveat 
pertaining to the Bayesian approach is that the final distribution of msusY depends on the prior 
distribution of msuSY in the most generic landscape mmm- As we have already mentioned, 
we assume that the prior distribution is not strongly biased towards high energies, although we 
do not need to make any particular assumptions regarding the prior distribution in the following 
argument. 

We should also mention that there have been several attempts in the literature to solve the 
question of the most probable SUSY scale by imposing further restrictions on the landscape 
of vacua for habitability. For example, the habitability condition has been used to restrict the 
landscape of vacua based on the abundance of dark matter in [ 22l [23l 23] . which leads to a 
sharp lower cut-off for the distribution of msuSY jf] In this paper, we shall refer to this type of 
restriction of the landscape based on the requirement of habitability as cosmological selection. 

Now, let us discuss how the domain wall problem related to i?-symmetry breaking provides 
us with a means of cosmological selection. In the above argument, we have eventually restricted 
the vacuum landscape so that all vacua in the landscape have the same electroweak scale uew- 
Before applying this restriction, let us now hypothesize that i?-symmetry breaking is caused by 
gaugino condensation in our vacuum. Then, too small values of the gravitino mass lead to unin¬ 
habitable universes for a given Hubble parameter during inflation (see Fig. [I]). Correspondingly, 
the distribution of msusY should be cut off for msuSY < m \/ 2 C^O according to cosmological 
selection for habitable universes (see Eq. ([ 3 ])) (see Fig.[3])(^] Finally, after applying the constraint 
on the electroweak scale, we obtain the distribution of mgusY in the landscape for a given uew- 
The crucial difference from the previous result is that the resultant distribution of msusY does 
not necessarily peak at uew anymore. Instead, it now peaks at my 2 (H) for my 2 (H) uew- 
In this case, the most probable sparticle masses can be much higher than uew> which gives us 
an answer why nature would prefer low-scale SUSY, but end up with a SUSY scale that is not 
completely natural. 

Interestingly enough, mgusY in the multi-TeV range, which is also suggested by the observed 
Higgs boson mass, can be reconciled with the concept of naturalness if the Hubble scale during 
6 See also [25] for a related discussion. 

‘Here, we have assumed that SUSY breaking is mediated to the MSSM sector via gravity mediation, i.e., 
msusY — m- 3 / 2 . 
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(a) Landscape of vacua 


(b) Landscape of vacua 
of the flat universe 


(c) Landscape of vacua with 
gaugino condensation 


(d) Landscape of vacua 
of v E w — 174GeV 



Vo ~0 


gaugirto v E w~174GeV 

condensation ■. • 1 • '. 


Distribution of msusY 


msusY 



Distribution of msusY 



Distribution of msusY 



vew« 


msusY 


Figure 3: Schematic illustration of how the probable range of tosusy can be deduced by means of 
cosmological selection under the assumption that the spontaneous breaking of R-symmetry is provided 
by gaugino condensation. From (a) to (b), we restrict the landscape to vacua with almost vanishing 
cosmological constant, which leads again to a distribution of tosusy sharply biased towards low energies. 
From (b) to (c), we restrict the landscape to vacua where the spontaneous breaking of R-symmetry is 
provided by gaugino condensation, which cuts off tosusy m.g, 9 (//). From (c) to (d), we restrict the 
landscape to vacua with an electroweak scale of vew — 174 GeV, which results in a distribution of tosusy 
that peaks around m^, 2 (H) ^w- 


inflation is in the range of 10 14 GeV (see Fig.[l|. As mentioned earlier, such a large Hubble 
parameter during inflation can be tested via future measurements of the tensor fraction in the 
CMB fluctuations. Therefore, the fundamental (and only seemingly metaphysical) question 
why nature would prefer low-scale SUSY, but end up with a SUSY scale that is not completely 
natural can be partially tested by future measurements of r. Conversely, if the tensor fraction is 
observed to be very small, the constraint on the gravitino mass from the domain wall argument 
becomes weak, and we will fail to reconcile the non-observation of sparticles with the concept 
of naturalness in this way. 

Several comments are in order. In the above argument, we have made the hypothesis that R- 
symmetry breaking is caused by gaugino condensation in our vacuum, which allows us to deduce 
the distribution of msusY hi consequence of cosmological selection. As we have mentioned in 
the previous section, there are, however, /7-symmetry breaking models which do not exhibit 
a domain wall problem. Thus, in order to fully answer the question whether the distribution 
of to-susy in a global landscape (for a given vew) really peaks at tosusy vew> we need to 
know the distribution of /7-breaking models, which goes far beyond the scope of this paper. 
The only thing we can say at this point is that we anyway need to live in a vacuum where 
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R-symmetry breaking is caused by gaugino condensation to reconcile the idea of naturalness 
with mguSY 3> vew by virtue of our domain wall argument. 

The same caveat applies to the distribution of the Hubble scale during inflation. That is, 
we need to know the distribution of H in a global landscape to conclude that the peak of the 
distribution of mguSY is really above uew- We only note here that chaotic inflation [26], which 
is free form the initial condition problem [27], predicts a Hubble scale of O(10 14 ) GeV during 
inflation^] The absence of the initial condition problem might explain why a large Hubble scale 
is chosen from the global landscape. Fortunately, our ignorance of the distribution of H can be 
compensated by future observations. That is, if our reasoning is correct, the tensor fraction of 
the CMB fluctuations will be measured to be rather large. 

In the above argument, we have shown that msuSY ^ fEW can be the most probable value 
when the Hubble parameter during inflation is high. However, we have not tried to explain why 
the electroweak scale is much smaller than msuSY- To answer this question, we expect that 
there are some anthropic reasons for the strength of the weak interaction, as in the case of the 
cosmological constant |30l l31j . In this paper, we, however, do not pursue these issues any further 
and. Instead, we refer to Refs. 133 133 EH EE] for more discussion on the anthropic selection for 
uew = 0(100) GeV. 

4 Summary and Discussion 

In this paper, we have discussed whether we can answer the fundamental question why nature 
would prefer low-scale SUSY, but end up with a SUSY scale that is not completely natural. This 
question is inevitable, if we postulate that low-energy SUSY is indeed realized in nature despite 
the null observation of sparticles at the LHC experiment below a TeV. This question becomes 
even more severe in view of the observed Higgs boson mass of about 125 GeV, which seems to 
point to sparticle masses in the multi-TeV range. As we have discussed, such a multi-TeV SUSY 
can be reconciled with the concept of naturalness under the assumption that the spontaneous 
breaking of an exact discrete R-symmetry is achieved by gaugino condensation in a pure SUSY 
Yang-Mills theory. In such theories, the dynamical scale of the Yang-Mills gauge interactions is 
required to be higher than the Hubble scale during inflation, in order to avoid the formation of 

8 For chaotic inflation models consistent with the constraint from the Planck satellite [2H|, see e.g. Refs. [25] 
and references therein. 
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domain walls, which puts a sharp lower cut on the distribution of the gravitino mass. With this 
sharp cut, we find that the distribution of tbsusy peaks in the multi-TeV range, if the Hubble 
parameter during inflation is of O(10 14 ) GeV. Our argument can be partially tested by future 
measurements of the tensor fraction in the CMB fluctuations. 

We should stress that what we have proposed in this letter is nothing less than a conceptual 
transition in how one should think about and address the big question of why SUSY has not yet 
been seen at the LHC. Conventionally, many people have attempted to construct models where 
the electroweak scale of 0(100) GeV is naturally obtained, even when mgusY lies in the multi- 
TeV. We, on the other hand, have taken a different approach in this letter, where we deduce the 
probable range of msxjSY for a given uew- We take the puzzle of the absence of sparticles at 
0(100) GeV as an important hint for the unknown structure of high energy physics. By adopting 
such a philosophy, we have, in fact, managed to infer the origin of R-symmetry breaking as well 
as the scale of the Hubble parameter during inflation in this paper. 

In the bulk of this paper, we have not made any assumption as to the sparticle spectrum of the 
MSSM. Let us comment here that our argument complies particularly well with a certain class 
of high-scale SUSY-breaking models where the gaugino masses are dominantly generated via 
anomaly mediation |36| (see also |37l 38] :39| for a discussion of the anomaly mediation mechanism 
in the superspace formalism of supergravity). In these models, no SUSY-breaking singlet fields 
are required and, hence, these models are free from the so-called Polonyi problem EDI SB 02]. 
Furthermore, the models in this class feature a good candidate for dark matter: the lightest 
gaugino (in particular the wino) in the TeV range or below. Therefore, this class of models has 
advantages in cosmology, which might enhance the probability of these models of being actually 
selected according to cosmological selection. Having gauginos in the TeV range (or below) is 
also important for the testability of the scenario. 

Throughout this paper, we have assumed an exact discrete R-symmetry. Inevitably, this 
symmetry should be anomaly-free mm- Related to this issue, let us consider the paradigm 
of pure gravity mediation (45j as an example. There, the R-charge of the MSSM Higgs bilinear 
H u Hd vanishes, and hence, a //-term of the order of the gravitino mass is naturally generated by 
the coupling of H u Hd to the VEV of the superpotential via Planck-suppressed operators E3S3 
(see also |48]L In this case, the difference between the MSSM contributions to the SU( 3) and the 
SU( 2) anomalies of the discrete R-symmetry is 4 and, hence, the exact R-symmetry is found 
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to be a Z 4 R symmetry]^] It should also be noted that an odd number of extra matter fields 
transforming in the 5 and 5* representations of SU(5) and with vanishing R-charge are required 
to make Z 4 R symmetry anomaly-free [35] , The existence of those extra matter fields, therefore, 
provides us with an additional possibility to test our assumption of an exact R-symmetry in 
future collider experiments. 

Finally, let us comment on the relation between R-symmetry and supersymmetric grand 
unified theories (GUTs). As shown in Ref. [52; [531 EH 155] , it is generically difficult to have an 
unbroken R-symmetry below the GUT scale in a class of GUT models where the standard model 


gauge groups are embedded in a single SU( 5) group 9 10 Thus, the existence of R-symmetry below 
the GUT scale fits well together with a class of GUT models where the standard model gauge 
groups are differently embedded into the subgroups of the GUT gauge group [56], 
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